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ABSTRACT 

We use observational data on the large scale structure (LSS) of the Universe mea- 
sured over a wide range of scales from sub-galactic up to horizon scale and on 
the cosmic microwave background anisotropics to determine cosmological parameters 
within the class of adiabatic inflationary models. We show that a mixed dark matter 
model with cosmological constant (AMDM model) and parameters flm = O.STIq'^si 



7+0.033 
-0.018' 



= 1.02 



+0.09 
-0.10' 



^A = o.eg+'azo' ^^ = 0.03+^:;;^, n, = i, n^ ^ o.os?^ 

h = 0.71^019, bci = 2.4lg 7 (Icr confidence limits) matches observational data on 
LSS, the nucleosynthesis constraint, direct measurements of Hubble constant, the high 
redshift supernova type la results and the recent measurements of the location and 
amplitude of the first acoustic peak in the CMB anisotropy power spectrum. The best 
model is A dominated (65% of the total energy density) and has slightly positive cur- 
vature, ri = 1.06. The clustered matter consists in 8% massive neutrinos, 10% baryons 
and 82% cold dark matter (CDM). The upper 2a limit on the neutrino content can 
be expressed in the form n^h? / N^-^'^ < 0.042 or, via the neutrino mass, nii, < 4.0eV. 
The upper l(2)cr limit for the contribution of a tensor mode to the COBE DMR data 
is T/S< 1(1.5). Furthermore, it is shown that the LSS observations together with the 
Boomerang (-f MAXIM A- 1) data on the first acoustic peak rule out zero- A models at 
more than 2a confidence limit. 

Key words: Cosmology: large scale structure - microwave background anisotropics 
- cosmological models: power spectrum - cosmological parameters 
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1 INTRODUCTION 

In the last decade of this century we have obtained impor- 
tant experimental results which play a crucial role for cos- 
mology: the COsmic Background Explorer has discovered 
the large scale ani sotropies of the cosm ic microwave back- 
ground radiation ( Bennett et al. 19961) ; the High-Z Super- 
nova Collaboration ( Riess et al. 19980 and the Supernova 
Cosmology Project (Perlmutter et al. 1998) found that the 
universe is accelerating rat her than decelerating ; the Super- 
Kamiokande experiment (Fukuda et al. 1998) discovered 



neutrino oscillations which prove the existence of neutri- 
nos with non-zero rest mass; balloon-borne measurements 
of the cosmic microwave ba ckground (CMB) temper ature 
fluctuations by Boomerang (le Bernardis et al. 2000 ) and 
MAXIMA- 1 ( [Hanany et al. 2000D have measured the height, 
position and width of the first acoustic peak which is in su- 
perb agreement with an adiabatic scenario of galaxy forma- 
tion. 

O n the other hand the comparison of recent experimen- 
tal dat.i on the large scale structure of the Universe with 



theoretical predictions of inflationary cosmology have shown 
since quite some time that the simplest cold dark matter 
(CDM) model is ruled out and we have to allow for a wider 
set of parameters to fit all observational data on the status 
and history of our Universe. These include spatial curvature 
(fifc), a cosmological constant (Oa), the Hubble parameter 
{h = iyo/(100km/s/Mpc)), the energy density of baryonic 
matter (^Ib), cold dark matter (flcdm), the number of species 
of massive neutrinos (A^^) and their density (fit/), the am- 
plitude of the power spectra of primordial perturbations in 
scalar (As) and tensor (At) modes and the corresponding 
power- law indices {us and nt), and the optical depth to 
early reionization (t). Constraining this multidimensional 
parameter space, determining the true values of fundamen- 
tal cosmological parameters, the nature and content of the 
matter which fills our Universe is an important and exciting 
problem of cosmology which has now become feasible due to 
the enormous progress in cosmological observations. About 
a dozen or more papers have been d evoted to this problem 
in th e last couple of years (see e.g. (Li neweaver & Barbosa 



199^ ; Lineweaver 1995 ; Efstathiou & Bond 199£ ; Pegmark 



© 0000 RA; 



2 R. Durrer, B. Novosyadlyj 



1999; [Bridle et al. 1999|; [Merchiorri et al. 2000|; [Novosyadly: 



E 
S 



et al. 2( [l00£J; [Tcgniark fc ZaldarriaRa 2000a|; rreRmark fc Zal 



darriaRJ 2000b|; [Lytli fc Covi 2000HNovosyadlyj et al. 200o'b[ 



Lange et al. 200C| ; [Balbi ct al. 2000t [Hu et al. 2000|; [TegmaTk 
et al. 2)001), some reviews are (pTirrcr fc Straumann 1999t 
Sahni fc Starobinsky 200C| ; [Priinack fc Gross 20Q(]| ; [Primack 
200C ) a nd references therein) . 



However, in spite of this intensive investigations the 
problem is still not satisfactorily resolved. Some of the re- 
maining issues are explained below. 

First of all, we would like to have observations which 
'measure' cosmological parameters in an as model indepen- 
dent way as possible. Clearly, most values of cosmological 
parameters obtained from observations of large scale struc- 
ture, galaxy clustering and CMB anisotropies are strongly 
model dependent. If the 'correct' model of structure forma- 
tion is not within the family investigated, we may not no- 
tice it, especially if the error bars are relatively large. This 
leads us to the next problem. Even if cosmological observa- 
tion have improved drastically, we still need more accurate 
data with better defined statistical properties (e.g we need 
to know the correlation of different measurements). The new 
CMB anisotropy data are already of this quality but the 
galaxy and cluster data are still relatively far from it. 

A next important point is the correspondence be- 
tween theoretical predictions and observational characteris- 
tics used in the analysis. We have to find a fast but accurate 
way to compute the theoretical values, especially when ex- 
ploring high dimensional parameter spaces. All parameters 
must be fitted simultaneously which renders the problem 
computationally complicated and very time consuming. Due 
to this difficulty, many authors search some subset of param- 
eters setting the others to some fixed 'reasonable' priors, 
thereby investigating a sub-class of cosmological models. As 
different authors also use different subsets of observational 
data, the resulting cosmological parameters still vary in a 
relatively wide range. 

Another problem are the degeneracies in parameter 
space which appear especial ly in the case when only C MB 
anisotropy data are used (Efstathiou & Bond 1999). It 



can be reduced substantially or even removed completely 
if galaxy clustering data, corresponding to different scales 
and redshifts, are combined with CMB measurements. This 
idea has already been employed on several occasions and is 
known under the name 'cosmic concordance' (for a recent 
review see Tegmark et al. (2000c)). 

The goal of this paper is to determine cosmological pa- 
rameters of the sub-class of models without tensor mode 
and no early reionization on the basis of LSS data related 
to different scales and diff'erent redshifts. In Novosyadlyj et 
al. (2000a) we have used the same approach to test flat 
models; we have shown that AMDM models are preferred 
in this class of models. There we have also shown that 
pure CDM models with h > 0.5, scale invariant primor- 
dial power spectrum, vanishing cosmological constant and 
spatial curvature are ruled out at very high confidence level, 
more than 99.99%. The corresponding class of mixed dark 
matter (M DM) models are ruled out at about 95% C.L. It 
was noted ( ISfovosyadlyj et al. 2000q ) that the galaxy clus- 
tering data set determines the amplitude of scalar fluctua- 
tions approximately at the same level as the COBE four- year 



data. This indicates that a possible tensor component in the 
COBE data cannot be very substantial. 

In this paper we test AMDM models with non-zero cur- 
vature. Furthermore, we use the data on the location and 
amplitude of the first acoustic peak determined from the 
most accurate recent measurements of the CMB power spec- 
trum. The data on the amplitude of the 2-nd and 3-d peaks 
is used as additional test for model prefered by large scale 
structure, COBE and 1-st peak data. We investigate the (in- 
)consistency of our data set with the 2-nd and 3-d peaks. We 
also use the SNIa constraint for comparison. 

The outline of the paper is as follows: in Sect. 2 we 
describe the experimental data set which is used here. The 
calculations of theoretical predictions and the method em- 
ployed to determine cosmological parameters are described 
in Sect. 3. In Sect. 4 we discuss our results and compare them 
with other investigations. Our conclusions are presented in 
Sect. 5. 



2 THE EXPERIMENTAL DATA SET 

Our approach is based on the quantitative comparison of 
the theoretical predictions for the characteristics of the large 
scale structure of the Universe with corresponding observa- 
tional ones. Theoretical predictions are calculated on the 
base of initial power spectrum of density perturbations of 
which shape strongly depends on all parameters supposed 
here for determination. So, model independent observational 
constraints on the inclination and amplitude of the power 
spectrum at different scales will be used in this search. 



2.1 CMB data 

We use the COBE 4-year dat a on CMB temperature 
anisotropies (Bennett et al. 1996) to normalize the density 
fluctuation power spectra according to Liddle et al. (1996) 
and Bunn & White (1997). Therefore, each model will match 
the COBE data by construction. 

We believe that using all available experimental data on 
AT/T at angular scales smaller than the COBE measure- 
ment is not an optimal way to search best-fit cosmological 
parameters due to their large dispersion (see for examples 
Fig. 10.1 of Durrer fc Straumann 1999, Fig. 2 of Novosyad- 
lyj et al. 2QQQa or Fig.l of Tegmark et al. 2000) which 
together with the large number of experimental points, 
~ 70 stipulating a high degrees of freedom, result in wide 
ranges for the co nfidence limits on cosmolo gical parameters. 
The Boomerang (|de Bernardis et al. 200C| ) and MAXIMA- 1 
( Hanany et al. 2000| ) experiments represent a new genera- 
tion of CMB measurements. They have produced a CMB 
map of about ~100deg^ with a resolution better than half 
a degree and a S/N~2, which allows to determine the loca- 
tion and amplitude of the first acoustic peak with high ac- 
curacy. The position of the first and amplitudes of the first, 
second and third acoustic peaks in the angular power spec- 
trum of the CMB temperature fluctuations together with 
the COBE data are the main measured characteristics of 
the CMB power spectrum. They contain information about 
amplitude and tilt of the primordial power spectrum of den- 
sity fluctuations at largest scales, from few tens of Mpc up 
to the current horizon scale of several thousand Mpc. They 
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are mai nly sensitive to the parameters Q,k, flmh^, ^a, i^th^ 



and to the normahzation of the initial power spectrum 



of density fluctuations. 

For example, the Boomerang data indicate that the first 
peak is located at the Legendre multipole £p = 197 ± 6 and 
has an amplitude of Ap = 69 ± 8/iK (this la error includes 
statistical and calibration errors). Here and in the following 
a tilde denotes observed quantities. We use these results in 
our search procedure. The MAXIMA-1 data {Ip « 220, Ap = 
78 ± 6/iK) marginally match Boomerang data and we will 
show that using them in combination with Boomerang data 
does not change the results significantly. The positions of 
2-nd and 3-d peaks are not well determined and we will 
not use them in the main sea rch procedure but we use their 
amplitudes as determined by ( Hu et al. 200C ) for comparison 
with the predictions of our best-fit model. 



2.2 Rich cluster data 

The important constraints on the form and amplitude of the 
matter power spectrum in the range from 10/i~^Mpc up to 
scales approaching 1000ft~'^Mpc can be obtained from the 
study of clusters of galaxies, their space distribution, mass 
and X-ray temperature functions. 

The power spectrum reconstructed from the observed 
space distribution of clusters has been determined many 
times for different samples from Abell, AGO and APM cat- 
alogs (see Einasto et al. 1997, Retzlaff et al. 1998, Tadros 
et al. 1998, Miller and Batuski 2000 and references therein). 
The remarkable feature of the determinations by different 
groups is similar slopes of cluster power spectra on scales 
0.02/iMpc"^ < fc < 0.1/iMpc"\ n ~ -1.5 (see above 
mentioned references). Here, we use the power spect rum of 
Abell- AGO clusters PA+ACo{kj) ( [Retzlaff et al. 1998| ) as ob- 
servational input. It is measured in the range 0.03/i/Mpc< 
k < 0.2/!,/Mpc where effects of nonlinear evolution are neg- 
ligible and it has well analyzed sources of uncertainties. 
The cluster power spectrum is biased with respect to the 
dark matter distribution. We assume that the bias is linear 
and scale independent. This is reasonable in the rang e of 
scales considered as predicted from local bias models (Fry 
& Gaztanaga 1995) and indicated by numeric al simula- 



tions (Benson ot al. 2000). In our previous paper (Novosyad 



lyj et al. 2QQQa) we have shown that not all the 13 points 
given in Retzlaff et al. (1998) are independent measurements 
and the effective number of degrees of freedom is 3. But to 
make best use of the observational information we use all 
13 points of the power spectrum to determine cosmological 
parameters and assign np = 3 for the number of degrees of 
freedom in the marginalization procedure. 

A constraint for the amplitude of the fluctuation power 
spectrum on cluster scale can be derived from the cluster 
mass and the X-ray temperature functions. It is usually for- 
mulated as constraint for the density fluctuation in a top-hat 
sphere of 8h~^ Mpc radius, as, which can be calculated for 
a given initial power spectrum P{k) by 

'^8 = TT^ / k^Pik)W^{8Mpck/h)dk, (1) 

2^^ Jo 

where H^(a;) = 3(sin j: — xcos a::)/^'' is the Fourier transform 
of a top-hat window function. Recent optical determinations 



of the m ass function of nearby galaxy clusters (Girardi et 
al. 1998|) give 



agfi^^ = 0.60 ± 0.04 



(2) 



where ai — 0.46 — 0.09Q,m for flat low-density models and 
Qi = 0.48-0.17^^ for open models (at the 90% G.L.). Sev- 
eral groups have found similar results using different meth- 
ods and different data sets (for a comprehensive list of ref- 
erences see Borgani et al. (1999)). This constraint on as is 
exponentially sensitive and thus allows only very small error 
bars. If the theory is correct this is of course a great advan- 
tage. However, if our understanding of cluster formation is 
not entirely correct, this will lead to discrepancies with other 
experimental constraints. 

From the observed evolution of the cluster X-ray tem- 
perature distribution function between z = 0.05 and z = 
0.32 we use the following constraint derived by Viana & Lid- 
die (1999) 



asflri 



0.56±0.19J7! 



0.11gnm+a2 



02 = 0.34 



for open models and 



(Tgf2° 



0.56±0.19J7! 



0.2 1gnm+Oi2 



a2 



0.47 



for flat models (both with 95% confldence limits). 

From the existence of three very massive clusters of 
galaxies observed so far at z > 0.5 an additio nal constraint 
has been established by (Bahcall & Fan 1998) 



as^° 



0.8 ±0.1 



(3) 



where as — 0.24 for open models and as — 0.29 for flat 
models. 

Note that all these constraints are given by slightly dif- 
ferent formulas for either JIa = or JIa + ^m = 1. However, 
we are going to use them for arbitrary values of JIa and 
Qm- Since our best fit models are relatively close to the fiat 
model, we mainly use the formula for the fiat case. We have 
checked that our results are insensitive to this choice. 



2.3 Peculiar velocity data 

Since our approach is based on the initial power spectrum of 
density fluctuations it seems most favorable to use the power 
spectrum reconstructed from the observed space distribu- 
tion of galaxies. But the galaxy power s pectra obtained from 
the two-dimensional APM survey (e.g. (Maddox et al. 199C; 
Tadros fc Estathiou 1996), and references therein), the GfA 

the 
and 



redshift survey (Vogclcy ct al. 1992; Park ct al. 199 



IRAS survey (Saunders ct al. 1992 



Saunders ct al. 200' 



the Las Gampanas Redshift Survey (ia Gosta et al. 1994 



Landy et al. 1996) differ signiflcantly in both, the ampli- 



tude and the position of the maximum. Moreover, nonlinear 
effects on small scales must be taken into account in their 
analysis. On the other hand, these power spectra contain 
large number of experimental points which are not indepen- 
dent and a decorrelation procedure for these power spectra 
must be employed. For these reasons and also in order to test 
the consistency between different data set, we do not include 
galaxy power spectra for the determination of parameters in 
this work. It will be interesting to compare our best flt pa- 
rameters with those obtained in analyzes including galaxy 
power spectra. 

Another constraint on the amplitude of the linear power 
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spectrum of density fluctuations in our vicinity comes from 
the stu dy of bulk flows of galaxies in spheres of large enough 



radii aiound our position. Since these data may be influ- 



enced by the local super-cluster (cosmic variance), we will 
use only the value of the bulk motion - the mean peculiar 
velocity of galaxies in a sphere of radius 50/i~^Mpc given by 
(kolatt & iSekel 19971), 



Fso = (375 ± 85)km/s. 



(4) 



With its generous error bars, this value is in a good agree- 
ment with other measurements of bulk mo tion at the scale 
40 - eOfe'^Mpc ( [Bertschinger et al. 1990 ; 



Courtoau et al 



1993 |; [qekel 1994|) (see also the review by Dekel 1999). 



2.4 Ly-a constraints 

An important constraint on the linear matter power spec- 
trum on small scales (fc ~ (2 — 40)/i/Mpc) comes from the 
Ly-a forest, the Ly-a absorption lines seen in quasar spec- 
tra (see Gnedin (1998), Croft et al. (1998) and references 
therein). Assuming that the Ly-a forest is formed by dis- 
crete clouds with a physical size close to the Jeans scale in 
the reionized inter-galactic medium at z ~ 2 — 4, Gnedin 
(1998) has derived a constraint on the value of the r.m.s. 
linear density fluctuations 



1.6 <aF{z^ 3) < 2.6 (95%C.L.) 
at kp ^ 34r2^/^/i/Mpc . 



(5) 



Taking into account the new data on quasar absorption lines, 
the effective equation of state and the temperature of the 
inter-g alactic medium at h igh redshift were re-estimated re- 
cently (Ricotti et al. 1999). As a result the value of Jeans 
scale at z; = 3 has moved to kp ~ S8Qm fe/Mpc (Gnedin 
200C). Here, we adopt this new value. 

The procedure to recover the linear power spectrum 
from the Ly-a forest has been elaborated by Croft et 
al.(1998). Analyzing the absorption lines in a sample of 19 
QSO spectra, they have obtained the following constraint 
on the amplitude and slope of the linear power spectrum at 
z = 2.5 and kp = 1.5Qrn h/Mpc, 



Ap(fcp) 



k;P{kp)/2n =0.57 ±0.26, 



(6) 
(7) 



at (Icr C.L.). The like constraints on the amplit ude and slope 
of the linear power spectrum was obtained by (McDonald et 



al. 2000) from the analysis of absorption lines in a sample of 



8 QSO. We will analyze these constraints in the context of 
our task and compare them with previous two. But in the 
main search procedure we will use the constraints given by 
Croft et al.(1998) as based on the more extensive sample of 
quasars. 

2.5 Other experimental constraints 

In addition to the CMB & LSS measurements described 
above we also use some results of global observations which 
are independent of the LSS model. For the value of the Hub- 
ble constant we set 



/i = 0.65 ±0.10, 



(8) 



which is a c ompromise between measur ements made by 



two grou ps, (rammann & Federspiel 1997) and (Madore et 
al. 1999). We also employ a nucleosynthesis constraint on 



the baryon density deduced from the determination of the 
primeval deuterium abundance 



Qbh'^ = 0.019 ± 0.0024 (95%C.L. 



(9) 



given by Buries et al. (1999) . The new, more precise deter- 
mination (Buries et al. 2000) confirms this value. 

Furthermore, we include the distance measurements of 
super novae of typ e la (SNIa) which constrain the cosmic 



expansion history ( Ries s et al. 1998; Perlmutter et al. 199J 



Perlmutter et al. 199S )). In a universe with cosmological 



constant this gives an important constraint on a combination 
of the values of the curvature, the cosmological constant and 
the matter content of the Universe. We use the following 



constraint in our parameter search ( Perlmutter et al. 1999| 



[n„ - 0.75!:!a] = -0.25 ± 0.125 



3 THE METHOD AND SOME TESTS 



(10) 



One of the main ingredients for the solution for our search 
problem is a reasonably fast and accurate determination of 
the linear transfer function for dark matter clustering which 
depends on the cosmological parameters. We use accurate 
analytical approximations of the MDM transfer function 
T{k; z) depending on the parameters flm, ^b, S7^, A^'^ and h 
by Eisenstein & Hu (1999). According to this work, the lin- 
ear power spectrum of matter density fluctuations is given 
by 



P{k; z) = Ask"^T'\k; z)DJ{z)/DJ{0), 



(11) 



where As is the normalization constant for scalar perturba- 
tions and Di{z) i s the linear growth factor, which can be 
approximated by (Carroll et al. 1992) 



Diiz) 



5a 

2~ 



.w 



l±z 



1 

70 



209Q^{z) - nljz) 4/7 



140 



+ ^Z'{z) 



where 

Q„,{Z) = n„(l + Z)V (f^m(l ± Zf + t^A ± f^fc(l ± Zf). 

We normalize the spectra to the 4-year COBE data 
which determine the amplitude of densitv perturbation a t 
horizon scale, 5h ( Liddle et al. 1996 ; Bunn and White 1997 ). 
The normalization constant As is then given by 



As = 27r^5^(3000Mpc//i)^+"= 



(12) 



The Abell-ACO power spectrum is related to the matter 
power spectrum at z = 0, P{k; 0), by the cluster biasing pa- 
rameter hci. As argued above, we assume scale-independent, 
linear bias 

PA+Aco{k)=hliP{k-Q). (13) 

For a given set of parameters ilm, Q.a, fit, Q.v, Ni,, Ua, 
h, and bd the theoretical values of PA+Acoik j) can now 



be obtain ed for the values kj (Table 1 of (Novosyadlyj et 
al. 2000a|)). We denote them by yj (j = 1, ...,13). 



The dependence of the position and amplitude of the 
first acoustic peak in the CMB power spectrum on cos molog- 
ical parameters has been investigated using CMBfast (Seljak 
fc Zaldarriaga 1996). As expected, and as we have shown in 
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our previous paper ( JNovosyadlyj et al. 2000a ), the results 
are, within reasonable accuracy, independent on the frac- 
tion of hot dark matter, fu = fi^/ilm, up to /^ 0.6. 

For the remaining parameters, Us, h, f2(,, flcd?n and JIa, 
we determine the resulting values £p and Ap using the an- 
alytical approximation given by Efstathiou & Bond (1999). 
We extend the approximation to models with non-zero cur- 
vature (Slfe = 1 — Q,n — SIa 7^ 0) by adding a coefficient for 
the amplitude and the peak location, which is determined 
numerically. The analytical approximation for the position 
of the first acoustic peak used here is 



tp = 0.7467r-v/3(l + Zr)—r-, :-, 

where tj* = Q^h^, and R = uj 

,,l/2sin(lc^fc|l/^y) 



(14) 



l/2 sinh(a;^^ y) 1/2 

Tj2 5 ^m Vi 



i-fc|i/= 



for open, flat and closed models respec- 



tively. Here y{LOm,i^k,i^A) is given by formula (8b) and 
Is{ujm,(jJb) by formulae (17-19) of Efstathiou & Bond 
(1999). The accuracy of this analytical approximation is bet- 
ter than 1%. 

Th a approaimation for the amplitude of first acoustic 



peak is LaaJaLkflaa 

At, = 



ep{ep + i) ^ r{ip + ^)rc^) 



2tt 



Tilp 



-|-0.838^(a;6 , ujcdm. , "-s )) 



1/2 



(15) 



where [nA{ujb,i^cd,n,ns) = 4.5(ns — 1) + ai + a-zUJcdm + 
asiOcdm + a-ii^l + asijJb + aaujbUJcdm + ajujk, with a\ = 2.376, 
a2 = 3.681, as = -5.408, 04 = -54.262, as. = 18.909, 
ae = 15.384, aj = 4.2 and C2 is the quadrupole anisotropy 
approximated by 



16 V c 



ns+3 



r(3 



V{2 + 



2 t 



p2(^ 



T{2 + - 



•-) 



(16) 



The values a\ — a% are the best-fit coefficients determined 
from a grid o f models computed with CMBfast (Efstathiou 



& Bone 1999). We have added the coefficient ay in order to 
account for curvature. The accuracy of Ap in the parameter 
ranges which we consider is better than 5%. We denote Ip 
and Ap by 3/14 and j/15 respectively. 

The theoretical values of the other experimental con- 
straints are obtained as follows: the density fluctuation ag 
is calculated according to Eq. (|l|) with P{k\ z) taken from 
Eq. (0. We set j/ie = as^m , yn = crs^m and j/is = (J&O.'S 
with corresponding values of a^ (i=l, 2, 3) for vanishing and 
non-zero curvature (see previous section). 

The r.m.s. peculiar velocity of galaxies in a sphere of 
radius R — 50/i~^Mpc is given by 



V50 = 



27r2 



fc^p(")me-'=^«/"'^ 



\k)e 



H/^(50Mpc k/h)dk, (17) 



where P^'"'{k) is power spectru m for the velocity fiel d 
of the density- weighted matter (Eisenstein & Hu 1999), 



VK(50Mpc k/h) is the top-hat window function. A previ- 
ous smoothing of the raw data with a Gaussian filter of 
radius Rf = 12/i~^Mpc is employed, similar to the proce- 
dure which has led to the observational value. For the scales 
of interest p("'(fe) Ri {0.°-^HofP{k;0)/k^. We denote the 
r.m.s. peculiar velocity by j/ig. 



The value by Gnedin (1998) from the formation of Ly- 
a clouds constrains the r.m.s. linear density perturbation at 
redshift z = 3 and wave number kp = SSQrn h/Mpc. In 
terms of the power spectrum, ap is given by 



^F(^)=2^y k'P{k; 



A-k/kpY 



dk. 



(18) 



It will be denoted by 2/20- The corresponding value of the 
constraint by Croft et al. (1998) is 



Al{kp,z) = klP{kp,z)/2n^ 



(19) 



at 2 = 2.5 with fcp = 0.008H{z)/{l + z)(km/s)~\ wiU be 
denoted by 2/21 ; 

H{z) = Ho [Qmil + zf + nfc(l + zf + Qa] ^'"^ 

is the Hubble parameter at redshift z. The slope of the power 

spectrum at this scale and redshift. 



Upiz) 



A log P(k.,z) 
A log k 



(20) 



is denoted by y-ii- 

For all tests except Gnedin's Ly-a clouds, we use the 
density weighte d transfer function Tcb,^{k,z) from (Eisen 



stein & Hu 1999). For Gnedin's ap we use Tcb{k, z) according 



to the prescription of Gnedin (1998). It must be noted that 
even in the model with maximal Q^ (~ 0.2) the difference 
between Tcb{k, z) and Tcbv(k, z) is less than 12% for k < kp. 
Early reionization changes somewhat the evolution of den- 
sity perturbation in the baryon component on small scales. 
This effect is not taken into account by the an alytical ap- 
proximation used here (Eisenstein & Hu 1999). Therefore, 



we restrict ourselves to models without early reionization. 
We calculate the Ly -a tests according to the prescription 
given in Sect. 5.4 of (Eisenstein & Hu 1999). 

Finally, the values of fit/i^, h and Qm — 0.75f2A are 
denoted by 3/23, 2/24 and 7/25 respectively. 

The squared differences between the theoretical and 
observational values divided by the observational error are 
given by x^, 



23 

E 



yj - yj 
Ay, 



(21) 



Here ijj and Aijj are the experimental data and their dis- 
persion, respectively. The set of parameters flm, ^a, ^v, 
Nu, 0,b, h, Us and bd are then determined by minimizing y 
using the Levenberg-Marquardt method (Press et al. 1992). 
The derivatives of the predicted values with respect to the 
search parameters which are required by this method are 
obtained numerically using a relative step size of 10~^ with 
respect to the given parameter. 

In order to test our method for stability, we have con- 
structed a mock sample of observational data. We start with 
a set of cosmological parameters and determine the "obser- 
vational" data for them which would be measured in case 
of faultless measurements with la errors comparable to the 
observational errors. We then insert random sets of starting 
parameters into the search program and try to recover the 
model which corresponds to the mock data. The method 
is stable if we can recover our input cosmological model 
(for more details of this test procedure see Novosyadlyj et 
al. (2000a). The code finds all the previously known parame- 
ters with high accuracy. Even starting very far away from the 
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true values, our method reveals as very stable and finds the 
'true' model whenever possible. This means that the code 
finds the global minimum of x independent of the initial 
values for the parameters. This also hints that our data set 
is sufficiently divers to be free of degeneracies (which plague 
parameter searches working with CMB data only). 



4 RESULTS AND DISCUSSION 



4.1 Calculations 



The determination of the parameters Qm, ^a, fi^, A^^, Qf,, 
h, Us an d hd by the Leven berg-Marquardt x^ minimization 
method (Press et al. 1992) can be realized in the following 
way: we vary the set of parameters ilm, ^a, ^v, ^b, h, Us 
and bci with fixed A^^ and find the minimum of x^ ■ Since 
A^i^ can be discrete we repeat this procedure three times for 
Nv = l, 2, and 3. The lowest of the three minima is the min- 
imum of x^ for the complete set of free parameters. Hence, 
we hav e seven free parameters. The formal number of ob- 



et al. 2000) and calculate them using the analytical approxi- 
mation given by the same authors then x^ « 18, which is far 
too much for 9 degrees of freedom. In this case the best-fit 
parameters are Q,^ = 0.37 ±0.07, ^a = 0.72 ±0.05, Q.^ ^ 0, 
Qt = 0.046±0.011, Us = 0.97±0.03, and h = 0.67±0.08. For 
the 2-nd acoustic peak and nucleosynthesis constraint the 
deviations of the predicted values from their observed coun- 
terparts are maximal {2.8a higher, and 1.4(7 higher respec- 
tively). If we exclude the nucleosynthesis constraint from 
the search procedure then x^ /Nf ~ 7/8 and best-fit param- 
eters become 0,^ = 0.34 ± 0.06, Q,a = 0.74 ± 0.05, Q.^ ^ 0, 
Qb = 0.055 ± 0.012, 71^ = 0.98 ± 0.03, and h = 0.72 ± 0.08. 
Practically all used constraints are satisfied but Qbh^ is 9(t 
higher then value deduced from th e determination o f the 
primeval deuterium abundance by ( Buries et al. 199!:) an d 
12o- higher then more recent value (Buries et al. 2000). 



servatic nal points is 25 but, as we have mentioned, the 13 
power spectra points can be described by just 3 degrees of 
freedom, so that the maximal number of truly independent 
measurements is 15. Therefore, the number of degrees of 
freedom for our search procedure is Np = A'oxp — A'par = 8 
if all observational points are used. In order to investigate 
to what extent the LSS constraints on fundamei ital param- 
eters matc h the constraints implied by SNIa (Perlmutter 
et al. 19991 ) we have determined all 8 parameters with and 
without the SNIa constraint (2/25). The results are presented 
in the Table |l[ 

Note, that for all models Xmm is in the range Nf — 
\/2Nf < Xmin ^ Nf + \/2Nf which is expected for a Gauss- 
ian distribution of Nf degrees of freedom. This means that 
the cosmological paradigm which has been assumed is in 
agreement with the data. (Note here, that the reduction of 
the 13 not independent data points of the cluster power spec- 
trum to three parameters is not important for our analysis 
since removing them from search procedure does not change 
the results essentially, as we will see later.) 

Let us investigate how the parameters of the best fit 
model vary if we include also the data of the MAXIMA- 1 
experiment. The location and amplitude of the first acous- 
tic peak determined fr om the combined Boomerang and 
MAXIMA-1 data are ( |Hu et al. 2000| ) (.j, = 206 ± 6, 



This problem of the inconsistency of the Boomerang and 
MAXIMA-1 values for the height of the second peak es- 
pecially with the nucleosynthesis constraint on the baryon 
abun dance has been dis c ussed at large in the recent litera- 
ture (Lange et al. 2000; Tegmark fc Zaldarriaga 2000b Hu 



et al. 2000t [Esposito et al. 2000| ; |Durrer et al. 200C| ). Since 



Ap = 78.6 ± 7. If we use them instead values used above, 
the best fit parameters remain practically unchanged, 0™, = 
0.37 ± 0.06, f^A = 0.66 ± 0.06, Q.^ = 0.03 ± 0.03, iV^ = 1, 
Q.b = 0.039 ± 0.010, Us = 1.05 ± 0.04, and h = 0.70 ± 0.09. 
Hence, including the MAXIMA-1 data into the determina- 
tion of the first acoustic peak is not essential in our anal- 
ysis and we will use here the values determined from the 
Boomerang data alone. This is however an important con- 
firmation of the consistency of the two data sets. 

We have also analyzed the influence of the amplitudes of 
the 2-nd and 3-d acoustic peaks on the determination of cos- 
mological parameters in the frame of our approach. If we add 
to our data set their values and errors as determined by (Hu 



we have nothing new to add to this subject here, we will 
not discuss it any further in this work. In what follows, we 
exclude the 2-nd and 3-d acoustic peaks from experimental 
data set in our search procedure but we will use them in the 
discussion of our best-flt model. 

The errors in the best-fit parameters as presented in 
Table hJ are the square roots of the diagonal elements of the 
covariance matrix which is calculated according to the pre- 
scription given in Press et al. 1992 (Chaper 15) or Tegmark 
& Zaldarriaga 2000a (Appendix A) . 



4.2 The best-flt model 

The model with one sort of massive neutrinos provides the 
best flt to the data, Xmin = 5.9. However, there is only a 
marginal difference in Xmin for A^i^ ~ 1, 2, 3. With the given 
accuracy of the data we cannot conclude whether massive 
neutrinos are present at all and if yes what number of degrees 
of freedom is favored. We summarize, that the considered ob- 
servational data on LSS of the Universe can be explained by 
a AMDM inflationary model with a scale invariant spectrum 
of scalar perturbations and a small positive curvature. 

Including of the SNIa constraint into the experimen- 
tal data set decreases Q.m, increases JIa slightly and prefers 
ill, ^ 0, a ACDM model. 

In Table we compare the values of the different ob- 
servational constraints with the predictions for the best-fit 
models (Table fy for N^ — 1). In both cases the calculated 
characteristics of the LSS are within the la error bars of 
the observed values. In the last row we indicate the age of 
the Universe determined according to the general e xpres- 
sion for non-zero cu rvature and non-zero A models (3ahni 
fc Starobinsky 200C|) 



^0 = Hq 



dz 



n™(i ± zY + f^fc(i + zY ± q,a{i + zY 



.(22) 



QW( 



'e treat Ha and Qm as free parameters, Q^ = 1 — Qa — ^n 



The predicted age of the Universe agrees well with recent 
determinations of the age of globular clusters. 

Comparing the results obtained without and with the 
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Table 1. Cosmological parameters determined from the LSS data described in the text without and with the SNIa constraint. The errors 
indicated are the square roots of the diagonal elements of the covariance matrix. 



N^ 



Ha 



Without SNIa constraint 



1 


5.90 


0.37±0.06 


0.69±0.07 


O.OSitO.OS 


0.037±0.009 


1.02±0.04 


0.71±0.09 


2 


6.02 


0.42±0.08 


0.64±0.09 


0.04±0.04 


0.038±0.010 


1.03±0.04 


0.71±0.09 


3 


6.17 


0.47±0.10 


0.59±0.08 


0.06±0.01 


O.OSSitO.OlO 


1.04±0.03 


0.70±0.09 



Including SNIa constraint 
0-3 6.02 0.32±0.05 0.75±0.06 < 10""' 0.038±0.010 1.0±0.05 0.70±0.09 



Table 2. Theoretical predictions for the used characteristics of the best-fit AMDM model with one sort of massive neutrinos with the 
cosmological parameters given in Tabic hi, first line (without SNIa constraint) and last line (including the SNIa constraint) are compared 
with observations. 





Observations") 




Predictions 




Characteristics 


Without SNIa 


Including SNIa 


^P 


197±6 


197 






197 


Ap 


69±8 


71.5 






71.9 


Vso, km/s 


375±85 


327 






308 


o-sf^m 


0.60±0.022 


0.61 






0.60 


OH^'S 


0.56±0.095 


0.58 






0.58 


<T8f1™^ 


0.8±0.1 


0.69 






0.71 


ap 


2.0 ±.3 


1.9 






1.9 


A2(fcp) 


0.57 ±0.26 


0.56 






0.59 


np(fcp) 


-2.25 ±0.2 


-2.20 






-2.20 


h 


0.65±0.10 


0.71 






0.70 


n^h^ 


0.019±0.0012 


0.019 






0.019 


Qm - 0.75Qa 


-0.25±0.125 


-0.14 






-0.25 


to, Gyrs'') 


13.2 ±3.0'=), 11.5 ±1.5'*' 


12.6 






13.5 



^' all errors are ±lcr, *>) is not used in the search procedure, '^'(Carretta ct al. 1999), '^^ ([Chaboycr ct al. 199^ ) 



SNIa constraint, we conclude that the values of the funda- 
mental cosmological parameters Q,m, f^A and Q.k determined 
by the observational characteristics of large scale structure 
match the SNIa test very well. This can be interpreted as 
independent support of the SNIa result in the framework 
of the standard cosmological paradigm. However, in order 
to elucidate how LSS data constraint cosmological param- 
eters, we analyze further only the model obtained without 
the SNIa constraint. 

The best fit values of cosmological parameters deter- 
mined by LSS characteristicsTJ are Q,^ = 0.37 ± 0.06, SIa = 
0.69 ± 0.07, Q.„ = 0.03 ± 0.03, iV^ = 1, f^t = 0.037 ± 0.009, 
Us = 1.02 ± 0.04, and h = 0.71 ± 0.09. The CDM density 
parameter is flcdm = 0.30 ± 0.10 and fi^ — —0.06 ± 0.13. 
The neutrino content, which is compatible with zero is very 



' We still include the direct measurement of h and the nucleosyn- 
thesis constraint in the analysis. 



badly determined (100% error). The obtained value should 
be interpreted as an upper limit to the neutrino contribu- 
tion. Below we will discuss this upper limit in more detail. 

The value of the Hubble constant is close to the result 
by Madore et al. (1999) and Mould et al.(2000), somewhat 
higher than the directly measured value given in Eq. (pi). The 
spectral index coincides with the prediction of the simplest 
inflationary scenario, it is close to unity. The neutrino matter 
density Qi, — 0.03 corresponds to a neutrino mass of rrii, = 
9'iQ^h^ ~ 1.4 eV but is compatible with within la. The 
estimated cluster bias parameter bd — 2.36 ± 0.25 fixes the 
amplitude of the Abell-ACO power spectrum (Fig. hi) . 

Recently, it has been shown ( Novosyadlyj 1999| ) that 
due to the large error bars, the position of the peak of P{k) 
at fc ~ 0.05h/Mpc does not infiuence the determination of 
cosmological parameters significantly. Only the slope of the 
power spectrum on scales smaller than the scale of the peak 
is relevant for cosmological parameters. On the other hand, 
the relation of the peak in PA+ACo{k) obtained from the 
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Table 3. Best-fit values of cosmological parameters determined from tiie different data set. 



No 



Data set 



^minl ^ 



Ua 



1 All observable data points are used 

2 PA+ACo{^y^ points are excluded 

3 £p, Ap are excluded 

4 V50 is excluded 

5 as^m^ is excluded 

6 ag^m^ is excluded 

7 (TsQrn^ is excluded 

8 All erg tests are excluded 

9 the first Ly-« test is excluded 

10 The second Ly-o test is excluded 

11 Both Ly-« tests are excluded 

12 The nucleosynthesis constraint is excluded 

12 The direct constraint on h is excluded 

13 Both previous constraints are excluded 

14 Vso, agQj^^ and Ap(fcp) are excluded 

15 5-8^^^ and Ap{kp) are excluded 



5.90/7 


0.37 


0.69 


0.027 


0.037 


1.02 


0.71 


2.12/4 


0.32 


0.75 


0.0 


0.039 


1.00 


0.70 


4.79/5 


0.39 


0.47 


0.058 


0.042 


1.14 


0.67 


5.54/6 


0.37 


0.69 


0.021 


0.038 


1.00 


0.71 


4.58/6 


0.45 


0.61 


0.052 


0.039 


1.03 


0.69 


5.88/6 


0.37 


0.69 


0.027 


0.037 


1.02 


0.71 


4.72/6 


0.38 


0.68 


0.028 


0.038 


1.01 


0.70 


3.85/4 


0.49 


0.57 


0.060 


0.041 


1.04 


0.68 


5.46/6 


0.42 


0.65 


0.048 


0.039 


1.02 


0.70 


5.81/5 


0.37 


0.69 


0.026 


0.037 


1.02 


0.71 


4.49/4 


0.56 


0.50 


0.21 


0.042 


1.04 


0.67 


4.52/6 


0.29 


0.89 


0.023 


0.001 


1.04 


0.67 


4.18/6 


0.29 


0.71 


0.038 


0.023 


1.05 


0.91 


4.16/5 


0.29 


0.71 


0.041 


0.013 


1.07 


0.87 


5.52/4 


0.37 


0.69 


0.021 


0.038 


1.01 


0.70 


5.88/5 


0.37 


0.68 


0.028 


0.037 


1.02 


0.71 
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Figure 1. The observed Abell-ACO power spectrum (filled cir- 
cles) and the theoretical spectra predicted by closed AMDM mod- 
els with parameters taken from Table n {N,^ = !)• 



Figure 2. The CMB power spectra nredicted by best-fit AMDM 
models with parameters from Table | {N^, = 1) and COBE DMR 
(Bennett et al. 1996), Boomerang (de Bernardis et al. 2000) and 
MAXIMA-1 (Hanany et al. 2000) experimental data. 



space distribution of Abell - AGO clusters around us to the 
matter density of the power spectrum of entire Universe is 
still under discussion. Using numerical simulations, Retzlaff 
et al. (1997) have shown that the pronounced peak in the 
spectrum (the fifth data point in Fig. hi) could be purely due 
to cosmic variance. Therefore, it should not influence cosmo- 
logical parameters. In fact, the maximum of our fitting curve 
is at a different position, which shows that this peak position 
is not relevant for the present work. The oscillation of the 
PA+ACo{k) around the best-fit P{k) in Fig. nl determined 
from all observable data on LSS reflects the real distribution 
of rich clusters of galaxies in the vicinity of ~ 300h~^Mpc 
of our own galaxy only. This is supported by similar fea- 
tures in spectra recons tructed from the expan ded sample 
of Abell-ACO clusters ([Miller fc Batuski 2000|) and IRAS 



Point Source Catalog Redshift Survey ( Saunders et al. 2000 



Hamihon et al. 2000) 



spectra of CMB temperature fluctuations for both best- 
fit models. Comparison with the COBE, Boomerang and 
MAXIMA-1 experiments are shown in Fig. ||. The CMB 
power spectrum predicted by both best-fit models matches 
the data very well within the range of the first acoustic 
peak. But it does not reproduce the absence of a second 
peak inferred from the Boomerang and MAXIMA-1 data at 
£ > 350. This problem has b een discussed intensively in lit- 
erature (Lange et al. 2000; Tegmark fc Zaldarriaga 2000b; 



Using CMBfast we have calculated the angular power 



Hu et al. 2000| ; [Esposito et al. 20001 ). The lack of power 
in this range strongly favors models with more baryons 
than compatible with standard cosmological nucleosynthe- 
sis. The MAXIMA-1 data reduces the problem so mewhat 
but does not remove it entirely ( Hu et al. 2000 ). How- 
ever, as we shall discuss, the cosmological parameters which 
match Boomerang and MAXIMA-1 data at high spheri- 
cal harmonics also strongly disagree with other LSS con- 
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straints used here (see Subsection 4.8 below). Furthermore, 
the Boomerang, MAXIMA- 1 and other CMB data in this 
range do not match each other very well. This (and the 
amount of work already published on this subject some of 
which is cited above) prompted us to ignore the problem of 
the second peak in the CMB anisotropy spectrum in this 
work. Future flights of Boomerang and MAXIMA and/or 
the future projects MAP and Planck will certainly decide 
on this very important issue, but we consider it premature 
to draw very strong conclusions at this point. 

Finally let us mention some global characteristics of a 
Universe with our best-fit cosmological parameters. Its age 
of io = 12.6 Gyrs is in the r ange of values determ i ned from 
the age of g lobular clusters ( Chaboyer et al. 1998 ; Carretta 
et al. 1J99). The deceleration parameter is qo = —0.52, in 
good agree ment with the SNIa con straint presented above 
leading to ( |Perlmutter et al. 1998| ) go = -0.57 ± 0.17. The 
original deceleration {q > 0) changes into acceleration {q < 
0) at the redshift Zd ~ 0.55. The 'equality epoch', pm{ze) = 
ph{ze), corresponds to the redshift Ze ~ 0.23. 



4.3 The influence of different experimental data 

One important question is how sensitive our result responds 
to each data point. To estimate this, we exclude some data 
points from the search routine and re-determine the best-fit 
parameters. The results of this procedure are presented in 
Table 0. In all cases when data on the first acoustic peak 
are included Q.m + ^a ~ 1.06, very slight positive curvature 
(fifc « —0.06) but compatible with flat, i.e. the geometry 
is defined mainly by the position of the flrst acoustic peak. 
The LSS data without CMB measurementsprefer an open 
Universe with Q.k = 0.14 (4 row in the Table ^. The value of 
fim never exceeds 0.56, Q,a is always larger 0.47 and in most 
cases Q.A > ^m- The best-flt values of the spectral index Us 
and h for the different observable data sets are in the rela- 
tively narrow ranges of, 0.99-1.14 and 0.67-0.72 respectively. 
The baryon content, Qb is fixed by the nucleosynthesis con- 
straint. Without this constraint (12 row in Table H) Qi, is 
lower, Qb ~ 0.001, even below the value inferred from the 
luminous matter in the Universe, f2ium ~ 7 x 10"'^. 

The hot dark matter content, Sli/, is reduced mainly by 
the Ly-Q constraints but it is poorly determined in all cases. 
If instead of or together with these Ly-a constraints we use 
those by McDonald et al. (2000) which reduce the power at 
small scales, then the best-flt value for the neutrino content 
is « 0.07. But in this case the predictions for Ly-a con- 
straints by Gnedin (1998) and Croft et al. (1998) are out of 
their la ranges. Moreover, the constraints by McDonald et 
al. (2000) are not in good agreement with other data, espe- 
cially, Bahcall & Fan (1998) and the SNIa constraints. We 
have not included these constraints any further in our deter- 
mination of cosmological parameters. Note however that the 
neutrino content is mainly constrained by the Ly-a data. If 
both Ly-a tests are excluded, the best fit value of fl^ raises 
to 0.21! 

Excluding the direct measurement of the Hubble pa- 
rameter from our search procedure leads to a substantially 
larger value of ft ~ 0.91 which is in disagreement with the 
direct determination. 

The comparison of the 1-st and 2-nd rows of Table H 
shows that the Abell-ACO power spectrum prefers a slope 



Table 4. The best fit values of all the parameters with errors 
obtained by maximizing the (Gaussian) 68% confidence contours 
over all other parameters. 



central value and errors 
parameter without SNIa constraint with SNIa constraint 



flm 


0.37 


f^A 


0.69 


n. 


0.03 


n^ 


0.037 


Us 


1.02 


h 


0.71 


bcl 


2.4- 



-HO. 25 

0.15 
-1-0.15 

0.20 
+0.07 

0.03 
+0.033 

0.018 
+0.09 

0.10 
+0.22 

0.19 
+0.7 

0.6 



0.32 
0.75 
0.0 
0.038 
1.00 
0.70 
2.2 



+0.20 

0.11 
+0.10 

0.19 
+0.09 

0.0 
+0.033 

0.019 
+0.13 

0.10 
+0.28 

0.18 
+0.8 

0.5 



of the matter power spectrum in the range 0.02 < fc < 
0.1/i/Mpc n ~ —1.5 which results in lowering JIa and in- 
troduces a small but non-zero neutrino content. 

The constraints agQ^^ (Viana & Liddle 1999) and 



Ap{kp) have practically no infiuence on the determination 
of parameters (rows 6, 10 and 15) due to their large error 
bars. They can be removed from the data set which reduces 
the number of effective degrees of freedom to Nf = 5; this 
is important for the marginalization procedure. 



4.4 Marginalization 

The next important question is: which is the confidence 
limit for each parameter marginalized over the others. The 
straight forward answer is the integral of the likelihood func- 
tion over the allowed range of all the other parameters. But 
for a 7-dimensional parameter space this is computationally 
time consuming. Therefore, we estimate the la confidence 
limits for all parameters in the following way. By varying all 
parameter we determine the 6-dimensional x^ hyper-surface 
which contains 68.3% of the total probability distribution. 
We then project this hyper-surface onto each axis in param- 
eter space. Its shadow on the parameter axes gives us the la 
confidence limits for the cosmological parameter under con- 
sideration. The la confidence limits obtained in this way for 
AMDM models with one sort of massive neutrinos are given 
in Table W. Including agQ^^ and Ap(fcp) does not change 
the marginalized limits significantly. 

It must be noted that even though the upper la edge 
for h is 0.93 when marginalized over all other parameters 
for the data used here, the resulting age of the Universe is 
still larger than the lowest value allowed for the age of the 
oldest globular clusters, to ~ 10 Gyrs if Qa > 0.72. In the 
Fig. we present the constraints in the Qa — h plane given 
by the lower limit for the age of the Universe, 10 Gyrs, for 
models with zero- and positive curvature. The range above 
corresponding line is excluded by this limit. Thus, the lower 
limit for the age of the Universe additionally constrains the 
confldence limits on the parameters, h and 0^ from above 
and on Qa from below. 

We have repeated the marginalization procedure includ- 
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Table 5. The upper limits for the neutrino content and mass (in 
eV) at different confidence levels. 





la 


C.L. 


2a 


C.L. 


3a 


C.L. 


N^ 


n. 


nil, 


n. 


niu 


n. 


nil, 


1 


0.10 


3.65 


0.13 


3.96 


0.18 


4.04 


2 


0.15 


2.79 


0.21 


3.06 


0.29 


3.35 


3 


0.20 


2.40 


0.27 


2.67 


0.35 


2.78 



0.4 0.6 

"a 



Figure 3. The lines in the Ha — h plane corresponding to the 
lower limit on age of the Universe of lOGyrs established from old- 
est globular cluster for models with zero- and positive curvature. 
The range below the corresponding line is allowed. 



ing the SNIa test (last column in Table 0). In this case we 
have to use all input data points (15 independent measure- 
ments), since neglecting asQ,"^ and Ap(fcp) does somewhat 
change the marginalized limit. Hence, the number of degrees 
of freedom is Np = 8 {la confidence limits corresponding 
to x^ ^ 15.3). The SNIa test reduces the confidence ranges 
of flm and Qa in spite of the larger number of degrees of 
freedom, but it results in somewhat wider la error bars for 
the other parameters due to the increase of Np and the low 
accuracy of the added data points. 



4.5 1 he status of some subclasses of models 



The errors shown in Table ^ define the range of each pa- 
rameter within which by adjusting the remaining parame- 
ters a value of Xmin < 11-8 can be achieved. Of course, the 
values of the remaining parameters always lay within their 
corresponding 68% likelihoods given in the Table W. Models 
with vanishing A are outside of this marginalized la range 
of the best-fit model determined by the LSS observational 
characteristics used here even without the SNIa constraint 
(column 2). Let us investigate the status of these models in 
more detail. For this, we set JIa = as fixed parameter and 
determine the remaining parameters in the usual way. The 
minimal value of x^ is X^ ~ 24 with the following values 
for the other parameters: f2,„ — 1.15, Sl^ = 0.22, N^, = 3, 
Qt = 0.087, Us = 0.95, h = 0.47, b^i = 3.7 (ag = 0.60)). 
This model is outside the 2a confidence contour of the best- 
fit model for N^ = 3 (Table |I| without SNIa test). The 
experimental data which disagrees most with A = is the 
data on the first acoustic peak. If we exclude it from the ex- 
perimental data set, Xmm ~ 5.8 for an open model with fol- 
lowing best-fit parameters: Qrn = 0.48, Qv — 0.12, A^'^ — 1, 
nt = 0.047, Us = 1.3, h = 0.64, b^i = 2.5 {ag = 0.82)). 
This model is inside the la confidence contour of the best- 
fit AMDM model obtained without data on the first acoustic 
peak (row 3 of Table H) . The reason for this behavior is clear: 
the position of the 'kink' in the matter power spectrum at 
large scales demands a 'shape parameter' F = Q,mh^ ~ 0.25 
which can be achieved either by choosing an open model 
or allowing for a cosmological constant. The position of the 



acoustic peak which demands an approximately fiat model 
then closes the first possibility. 

Results change only slightly if instead of the Boomerang 
data we use Boomerang-fMAXIMA-l as discussed in Sec- 
tion 4.1. Hence, we can conclude that the LSS observational 
characteristics together with the Boomerang (-I-MAXIMA- 
1) data on the first acoustic peak already rule out zero-A 
models at more than 95% C.L. and actually demand a cos- 
mological constant in the same bulk part as direct measure- 
ments. We consider this a non-trivial consistency check! 

Flat A models in contrary, are inside the la contour 
of our best-fit model. Actually, the best fit fiat model has 
Xmin ~ 8.3 and the best fit parameters dm ~ 0.35 ± 0.05, 
^A = 0.65 T 0.05, n^ = 0.04 ± 0.02, A^^ = 1, Qt = 0.029 ± 
0.005, Us = 1.04±0.06, h = 0.81 ±0.06, b^i = 2.2±0.2 (gg ^ 
0.96) are close to our previous (Novosyadlyj et al. 2000a) 



results with a somewhat different observational data set. 

It is obvious, that fiat zero-A CDM and MDM mod- 
els are ruled out by the present experimental data set 
at even higher confidence limit than by data without the 



Boomerang and MAXIMA-1 measurements in (Novosyadlyj 
ct al. 2000 j ). 



4.6 Upper limits for the neutrino mass 

Since the neutrino content is compatible with zero, we deter- 
mine an upper limit for it. We first determine the marginal- 
ized la, 2a and 3a upper limits for fi^ for different values of 
N,^. Using the best-fit value for h at given Q^, we can then 
determine the corresponding upper limit for the neutrino 
mass, m„ = 94Qvh^ /Ni,. The results are presented in Ta- 
ble 0. For more species of massive neutrino the upper limit 
for f2^ is somewhat higher but rrii, is still lower for each C.L. 
The upper limit for fli, raises with the confidence level as 
expected. But the upper limit for the mass grows only very 
little due to the reduction of the best-fit value for h. The up- 
per limit for the combination Q,„h^ /N^'^'^ is approximately 
constant for all number species and confidence levels. The 
observational data set used here establishes an upper limit 
for the massive neutrino content of the universe which can 
be expressed in the form Q.^h^ /N°-'^'^ < 0.042 at 2a confi- 
dence level. The corresponding upper limit on the neutrino 
mass rriy < 4eV is close to the value obtained by (Croft et 
al. 1999|). 



4.7 Limiting the tensor mode 

Up to this point we ignored uncertainties in the COBE nor- 
malization. The statistical uncertai nty of the fit to the four - 
year COBE data, Sh, is 7% {la) (|Bunn and White 1997; 



© 0000 RAS, MNRAS 000, 000-000 



Cosmological parameters from complementary observations of the Universe 11 



and we want to study how this uncertainty influences the 
accuracy of cosmological parameters which we determine. 

Varying 5h in the la range we found that the best-fit 
values of all parameters except Q^, do not vary by more 
than 7% from the values presented in Table 1. Only f2^, on 
which Ict errors are of the order of 100%, varies in a range of 
20% . These uncertainties are significantly smaller than the 
standard errors given in Table 1 and ignoring them is thus 
justified. (Including this error raises our standard Ict errors 
from typically 10% - 20% to 11% - 21%.) 

Our results depend on a possible tensor component only 
via the COBE data which enters our calculation through 
the normalization constant Sh, in Eqns. (hlUld). We can es- 
timate the maximal contribution of a tensor mode in the 
COBE AT/T data in the following way: we disregard the 
COBE normalization and consider Sh as free parameter to be 
determined like the others. Its best-fit value then becomes 
6h^^ = (2.95 ± 2.55) • 10"^ (for iV^ = 1), while the best- 
fit values of the other parameters are Qm ~ 0.40 ± 0.08, 
^A = 0.66 ± 0.07, Q^ = 0.05 ± 0.05, fib = 0.038 ± 0.010, 
Us = 1.14 ± 0.31, h = 0.71 ± 0.09 and b^i = 2.4 ± 0.3. 
The best-fit value for density perturbation at horizon scale 
from the 4-year COBE data for this set of parameters is 
larger then the best-fit value determined from LSS char- 
acteristics, (5jr°^^ = 4.0 • 10"^ > 5h^^. This means that 
COBE AT/T data may contain a non-negligible tensor con- 
tribution. The most likely value of its fraction is given by 
T/S = (5^°^^ - (5^^^)/5^^s. This value is T/S = 0.36 for 
the corresponding best-fit values of 5^ a-nd 5^ from 

the Boomerang data alone and T/S = 0.18 from the com- 
bined Boomerang -f MAXIMA- 1 data. Since the standard 
error is rather large, ~ 90%, we determine upper confidence 
limits for T/S by marginalizing 5^ over all the other pa- 
rameters like we did for the neutrino content (see subsection 
4.6). We then obtain T/S < 1 at la C.L. and T/S < 1.5 at 
2ct C.L. from the Boomerang data alone for the amplitude 
and position of the first acoustic peak. If we use the com- 
bined Boomerang -f MAXIMA- 1 data these limits are some- 
what lower, 0.9 and 1.3 correspondingly, due to the higher 
amplitude of the first acoustic peak measured by MAXIMA- 
1. The Ict upper constraint on the tensor mode obtained 
recently by Kinney et al. (2000) from the Boomerang and 
MAXIMA- 1 data on the CMB power spectrum for the same 
class of models (T/S< 0.8 in our definition) is very close to 
the value obtained here. 



4.8 CHuiiipai'iyuii with utliyr parametbi' bstlmatluiiy 

The cosmological parameters determined here from 
LSS-I-CMB data agree well with the values obtained by 
other methods (see e.g. the review by Primack (2000)). The 
marginalized Ict ranges are still rather large due to the large 
experimental errors, the large number of parameters and the 
high degree of freedom. But this does, of course, not mean 
that an arbitrary set of parameters within the marginalized 
ranges matches the experimental data set with an accuracy 
< Ict. 

We compare our best fit model with others found in 
the recent literature by testing our data set as well as the 
Boomerang and MAXIMA- 1 data on the CMB power spec- 
trum. At first we calculate the predictions of the follow- 
ing models for our data set {Q,m,^A,^b,ns,h) = P = 



(0.49, 0.56, 0.054, 0.92, 0.65) obtained by Lange et al. (2000) 
as best-fit model for the Boomerang and LSS data (denoted 
there as model P9); P = (0.68, 0.23, 0.07, 1, 0.6) obtained by 
Balbi et al. (2000) as best-fit model to the MAXIMA- 1 and 
COBE DMR data; P = (0.35, 0.65, 0.036, 0.95, 0.8) obtained 
by Hu et al. (2000) as best-fit model to the Boomerang 
+ MAXIMA- 1 data on the first, second and third acoustic 
peaks; P = (0.3,0.7,0.045,0.975,0.82) obtained by Jaffe et 
al. (2000) as best-fit model to the Boomerang -I- MAXIMA- 
1 + COBE data on the CMB power spectrum; and the 
"concordance" model by Tegmark et al. (2000) which fa- 
vors P = (0.38,0.62,0.043,0.91,0.63). Some authors give 
several sets of parameters obtained for different priors or by 
including different data sets, we take the one from which 
we obtain a minimal x^ for our data set. All these models 
have no massive neutrino component, no tensor mode and 
reionization is either not included or can be neglected. The 
predictions of cosmologies with the above parameters for the 
data considered in this work are presented in Table pi 

The x^ presented in last row includes also Xj^+aco ~ 

L-li = l <y APAJrACOyki) j 

cluster bias, feci, which is considered as free parameters in 
each model. In spite of the fact that all parameters of each 
model are within the marginalized Ict ranges of the parame- 
ters of our best-fit model, the total value of x^ for the entire 
parameter sets rules out all the models at more than 2ct 
confidence level. Table indicates the crucial tests. Models 
A and C are ruled out mainly by the nucleosynthesis con- 
straint and the first cts test (cluster mass function). Model 
B strongly disagrees with all CTg tests (14ct, 2.6ct and 1.6ct 
correspondingly), both Ly-a tests (2.6ct and 2.5ct), the nu- 
cleosynthesis constraint (5.2ct) and the data on the location 
of the first acoustic peak (5.6ct). Moreover, models A and 
B do not match the SNIa test which we have not included 
into x^ . Model D strongly disagree with nucleosynthesis con- 
straint (9.4ct) and the Boomerang data on the location of the 
first acoustic peak (2.6ct). Model E does not match first and 
third CT8 tests (at 5. Ict and 2.4ct respectively), the first Ly- 
a test (at 2.4ct) and the data on the location of the first 
acoustic peak (3.2ct). The latter is due to the fact that the 
MAXIMA- 1 peak position is more than Ict away from the 
peak position derived by the Boomerang data alone. 

We now calculate the CMB power spectra for these 
models using CMBfast (version 3.2) and co mpare them with 
the expe rimental data from Boomerang (de Bernardis et 



al. 2000) and MAXIMA- 1 (Hanany et al. 2000). The x 



deviations for all models including our best-fit model are 
presented in the Table M. The first number indicates the x'^ 
for the range of the first acoustic peak, 50 < ^ < 375 (7 and 
5 data points for the Boomerang and MAXIMA- 1 exper- 
iments respectively), for the second number we have used 
the entire range, 50 < ^ < 750 (12 and 10 data points for 
the Boomerang and MAXIMA- 1 experiments respectively). 
In the range of the first acoustic peak our model fits as well 
as the other models, but the observed power spectrum at 
higher spherical harmonics is not reproduced by our model 
as we mentioned above. 

Therefore, models which match the Boomerang and/or 
MAXIMA- 1 CMB power spectrum at high spherical har- 
monics (in the range of second and third acoustic peak) dis- 
agree with some of the CTg, Ly-a and/or the nucleosynthesis 
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Table 6. Theoretical predictions for the observational values by best-fit models from the literature: A (Lange et al. 2000), B (Balbi et 
al. 2000), C (Hu et al. 2000), D (Jaffe et al. 2000), E (Tegmark et al. 2000). 





Observations 




Predictions 






Characteristics 


A 


B 


C 


D 


E 


^ 


197±6 


206 


231 


206 


213 


225 


A, 


69±8 


57 


68 


63 


72 


62 


V50, km/s 


375±85 


280 


310 


303 


293 


239 


CTSH^I 


0.60±0.022 


0.64 


0.91 


0.68 


0.58 


0.43 


'Tsn^^ 


0.56±0.095 


0.62 


0.93 


0.65 


0.65 


0.42 


CT8^m 


O.SitO.l 


0.70 


0.96 


0.79 


0.73 


0.49 


0"F 


2.0 ±.3 


1.7 


2.8 


2.2 


1.9 


1.1 


A2(fcj,) 


0.57 ±0.26 


0.51 


1.21 


0.81 


0.62 


0.25 


TlpyKpj 


-2.25 ±0.2 


-2.25 


-2.15 


-2.21 


-2.22 


-2.30 


h 


0.65±0.10 


0.65 


0.60 


0.80 


0.82 


0.63 


Qth^ 


0.019±0.0012 


0.023 


0.025 


0.023 


0.030 


0.02 


n,n - O.TS^A 


-0.25±0.125 


0.07 


0.51 


-0.13 


-0.23 


-0.09 


x' 




27 


285 


39 


105 


106 



Table 7. The x^ deviation of theoretical predictions for the CMB power spectrum from experimental results for the models in Table a 
and for our best-fit model. The first number represents the the value of x^ for the CMB power spectrum in the range of the first acoustic 
peak, 50 < £ < 375, the second number is for the entire range 50 < £ < 750. Clearly our model parameters are in serious disagreement 
with the experimental CMB data beyond the first acoustic peak. 



Experiment 



our best-fit model 



Boomerang 7.3/12.6 77.2/96.7 6.1/12.8 18.3/24.5 13.6/24.5 11.3/108.5 

MAXIMA-1 16.9/18.7 4.6/11.4 15.2/17.0 10.3/11.7 16.2/21.6 11.1/48.5 



constraints. And vice versa, model which match very well the 
LSS observational characteristics predicts the CMB power 
spectrum which disagrees with measurements by Boomerang 
and MAXIMA-1 on very small scales. The resolution of this 
problem can go in several directions. If the Boomerang and 
MAXIMA-1 measurements are confirmed, nucleosynthesis 
may have been more compl icated than assumed for the con- 
straint used in this work (Esposito et al. 2000). An other 



problem may be the cluster mass function constraint which 
is exponentially sensitive to the value of ag and might be 
too constraining, especially in view of all the uncertainties 
in the theory of cluster formation. Therefore, our constraint 
as^m = 0.60±0.022 has to be taken with a grain of salt and 
its incompatibility with, e.g. the CMB data may also hint to 
a problem in the theory of cluster formation. Last but not 
least, if inconsistencies in the determination of cosmological 
parameters persist even after a serious improvement of data, 
e.g. with the Sloan digital sky survey, this may hint that the 
correct model is not within the class considered. If we want 
to fit a snail within the class of all known mammals by x^ 
minimization (or by a much more sophisticated method) , we 
never obtain a very convincing fit. 



5 CONCLUSIONS 

The main observational characteristics on LSS together with 
recent data on the amplitude and location of the first acous- 
tic peak in the CMB power spectrum, and the amplitude of 
the primordial power spectrum inferred by the COBE DMR 
four year data prefer a AMDM model with the following pa- 

+0.07 

-0.03' 

+0.22 



rameters: Qm = 0.37+°i5, ^a = 0.69+°;2o, ^>. = 0.03 



bet 



- 1, fti 
2 4+°-^ 



Q QO7 + O.O33 
U.UO/_gQ^g, rig 



1 02+°''^ h 



0.71 



1(7 marginalized ranges). 

The central values correspond to a slightly closed {Qk ~ 
—0.06) AMDM model with one sort of 1.4eV neutrinos. 
These neutrinos make up about 8% of the clustered mat- 
ter, baryons are 10% and the rest (82%) is in a cold dark 
matter component. The energy density of clustered matter 
corresponds to only 35% of the total energy density of mat- 
ter plus vacuum which amounts to f2 = 1.06. The massive 
neutrino content is compatible with zero and we have es- 
tablished an upper limit in the form of n^h'^/N°-^^ < 0.042 
at 2a confidence level. The upper 2a limit for the neutrino 
mass is 4.0eV. 

If COBE normalization is disregarded, the best-fit value 
of the density perturbation at horizon scale is Sf^ = (2.95± 
2.55) • 10~^ while the best-fit values of the other parameters 



© 0000 RAS, MNRAS 000, 000-000 



Cosmological parameters from complementary observations of the Universe 13 



are Q.^ = 0.40, ^a = 0.66, ^^ = 0.05, iV^ = 1, fib = 0.038, 
Us — 1.14, h — 0.71 and bd — 2.4. Comparison it with 
the best-fit value to the COBE 4-year data 5^ gives an 

estimate for the contribution of a tensor mode to the COBE 
DMR data: T/S = 0.36to;|6 from the Boomerang data on 
the first acoustic pealc and T/S = O.lSlgjg {la confidence 
hmits) when the combined Boomerang-l-MAXIMA-l data 
are used. The upper limits on T/S at 2a C.L. for these two 
cases are 1.5 and 1.3 respectively. 

The values for the matter density Qm and the cosmo- 
logical constant JIa for the best-fit model are close to those 
deduced from the SNIa test. Including this test in the ob- 
servational data set, results to a somewhat larger value of 
f^A (7%) and slightly lowers Qm- 

The observational characteristics of large scale structure 
together with the Boomerang (+MAXIMA-1) data on the 
first acoustic peak rule out zero-A models at more than 2a 
confidence limit. 
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